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Abstract 

Thermodynamic properties of the QCD nonperturbative vacuum with two hght quarks 
are studied. It is shown that at low temperatures T < relativistic massive pions can be 
treated within the dilute gas approximation. Analytic temperature dependence of the quark 
condensate is found in perfect agreement with the numerical calculations obtained at the 
three-loop level of the chiral perturbation theory with non-zero quark mass. The gluon con- 
densate slightly varies with the increase of the temperature. It is shown that the temperature 
derivatives of the anomalous and normal (quark massive term) contributions to the trace of 
the energy-momentum tensor in QCD are equal to each other in the low temperature region. 
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1. The investigation of the vacuum state behavior under the influence of various external factors 
is known to be one of the central problems of quantum field theory. In the realm of strong interac- 
tions (QCD) the main factors are the temperature and the baryon density. At low temperatures, 
T < Tc { Tc-temperature of the " hadron-quark-gluon" phase transition ), the dynamics of QCD is 
essentially nonperturbative and is characterized by confinement and spontaneous breaking of chiral 
symmetry (SBCS). In the hadronic phase the partition function of the system is dominated by the 
contribution of the lightest particles in the physical spectrum. In QCD this role is played by the 
TT-meson which is the Goldstone excitation mode in chiral condensate in the limit of two massless 
flavors. Therefore the low temperature physics can be described in terms of the effective chiral 
theory jl], ^, ^ . An important problem is the behavior of the order parameter (the quark conden- 
sate (qq)) with the increase of the temperature. In the ideal gas approximation the contribution 
of the massless pions into the quark condensate is proportional to |§, § • In chiral perturbation 
theory (ChPT) the two - and three-loop contributions and ~ correspondingly) into (qq) 

have been found in |^ and Q.In the general case of Nf massless flavors, the low temperature 
expansion of the quark condensate takes the form 



A„ 



~J^-' ^[Y^J "''^^''f''\l2Fl) + W 

where ~ 93MeV is the value of the pion decay constant and A, ~ 470Mev {Nf = 2)0. The 
result (0) is an exact theorem of QCD in chiral limit. 

The situation with the gluon condensate (G^) = ((gGJ^j,)^)is very different. The gluon conden- 
sate is not an order parameter of the phase transition and it does not lead to any spontaneous 
symmetry breaking (SSB). At the quantum level the trace anomaly leads to the breaking of the 
scale invariance which is phenomenologically described by non-zero value of (G^). However, this 
is not a SSB phenomenon and hence does not lead to the appearance of the Goldstone parti- 
cle. The mass of the lowest excitation (dilaton)is directly connected to the gluon condensate, 
mo cx ((G^))^/"'. Thus, in the gluonic sector of QCD, the thermal excitations of glueballs are 
exponentially suppressed by the Boltzmann factor ^ exp{— TOgj/T} and their contribution to the 
shift of the gluon condensate is small {A{G^)/{G^) -0.1 % at T -200 MeV) |||. Next we note that 
in the one-loop approximation of ChPT pions are described as a gas of massless noninteracting 
particles. Such a system is obviously scale-invariant and therefore does not contribute to the trace 
of the energy-momentum tensor and correspondingly to (G^). As it has been demonstrated in joj 
the gluon condensate temperature dependence arises only at the ChPT three-loop level due to the 
interaction among the Goldstone bosons. 

For the quark condensate the account of the non-zero quark mass was done in 0. Within ChPT 
two masses of the u and d quarks are considered as perturbation (m„ and rrid are much smaller 
than any scale of the theory). Correspondingly thermodynamic quantities ( pressure, {qq){T)) are 
represented as a series in powers of T and m in low temperature expansion. As a consequence of 
the fact that the spectrum of the unperturbed system contains massless pions, the quark mass term 
generates infrared singularities involving negative powers of T . In Ref . [Q results of a numerical 
procedure was presented in the form of the graphical dependence of (qq) {T) as a function of T, 
which we will compare with our analytic formula. The deviation of the quark condensate from the 
chiral limit (]l|) turns out to be substantial (~ 25%) at T = 140 MeV in both cases. 

It is well known that due to the smallness of pion mass as compared to the typical scale of strong 
interactions, the pion plays a special role among other strongly-interacting particles. Therefore 
for many problems of QCD at zero temperature the chiral limit, Af^ — > 0, is an appropriate one. 
On the other hand a new mass scale emerges in the physics of QCD phase transitions, namely the 
critical transition temperature Tc- Numerically the critical transition temperature turns out to be 
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close to the pion mass, Tc ~ M^rj^ Thus m the confining phase of QCD the smaU parameter M^^/T 
is lacking and hence the interval <^ T < Tc where the high temperature expansion is applicable 
for thermal pions does not exist. However hadron states heavier than pion have masses several 
times larger than Tc and therefore their contribution to the thermodynamic quantities is damped 
by Bolzmann factor ~ eicpl— Mhadr /T} . Thus the thermodynamics of the low temperature hadron 
phase, T < Af^, is described basically in terms of the thermal excitations of relativistic massive 
pions. 

In the present paper we study the behavior of the quark and gluon condensates in the low 
temperature phase. It is shown that the relativistic dilute massive pion gas approximation is 
applicable at temperature T < M-^. Analytical expressions for condensates are obtained. The 
comparison with the results of three-loop ChPT calculations {qq){T) and {G'^){T) is made. The 
low temperature relation for the trace of the energy-momentum tensor in QCD with two light 
quarks is obtained based on the general dimensional and renormalization-group properties of the 
QCD partition function and dominating role of the pion thermal exitations in the hadronic phase. 

2. For non-zero quark mass (m,q ^ 0) the scale invariance is broken already at the classical 
level. Therefore the pion thermal excitations would change, even in the ideal gas approximation, 
the value of the gluon condensate with increasing temperature. To determine this dependence can 
be use the general renormalization and scale properties of the QCD partition function. This is 
a standard method and it is used for derivation of low-energy QCD theorems For QCD at 
finite temperature and chemical potential these theorems were derived in [ p^ |13||. This method 
was used for investigation of QCD vacuum phase structure in a magnetic field |14| and at finite 
temperature and also in nuclear matter [ p^ . 

A relation between the trace anomaly and thermodynamic pressure in pure-glue QCD was 
obtained in [ p7[ by making use of the dimensional regularization in the framework of the renor- 
malization group (RG) method. Also within the RG method, but by employing slightly different 
techniques an analogous relation was derived in the theory with quarks in ref. p8[ . Then, by 
making use of the results obtained in the above-mentioned papers, we can write the following 
expression for the trace anomaly in QCD with quarks: 

where Pr is the renormalized pressure, /3(as) = das{M)/d InM is the Gell-Mann-Low function 
and 7m is the anomalous dimension of the quark mass. It is convenient to choose such a large 
scale that one can take the lowest order expressions, /3(as) — > —ba1/2Ti, where h = {llNc — 2Nf)/3 
and 1 + jrri 1- Thus, we have the following equations for condensates 

(0,r).?f!(4-rA-i:™,^)p„.i>p„. (3) 

3. The density of the pionic gas at temperature T is given by 

f d^p 1 
= V (27r)3 exp( VFTW/r) - r 

— 1/3 

and the average distance between particles in this gas is X^r — f^ir • The dilute gas approximation 
is valid, when the average distance L-j^ is much smaller than the mean free path Atttt calculated 

^The deconfining phase transition temperature is the one obtained in lattice calculations Tc(Nf = 2) ~ 173 MeV 
and Tc{Nf = 3) ~ 154 MeV O. 
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taking into account collisions, L^r ^ ^Ti-n- For a hadronic gas, the mean free path is well known 



[ [19| , 20 at low temperatures, where the gas almost exclusively consists of pions. An analysis of 
the TTTT collision rate shows that in the range 50 MeV< T <140 MeV, ChPT formula for mean 
free path A^^ ~ 12F^/T^ Q is vahd to within about 20% (2l). In Fig.l the behavior of the ratio 
L^/Atttt as a function of T is shown. Hence at T < the gas approximation for pions can be 
used. 




Figure 1: The ratio of the mean interparticle distance L-^ ( in pionic gas ) to the mean free path 
A,r7r as a function of the temperature 

Thus the effective pressure from which one can extract the condensates {qq){T) and {G'^){T) 
using the relations (||) and (Q) has the form 

Peff{T)^-e,ac+PAT). (6) 

where e^ac = is the nonperturbative vacuum energy density at T = and 

(M = -^(G'>+ E ^'i^m) (7) 

q—u,d 

is the trace of the energy-momentum tensor. In Eq.(^ P^r is the massive pions pressure 

= -^^Jwf ^""^^ ~ e'VP^JT^. (8) 
The quark and gluon condensates are given by the equations 

{G'){T} = DP.„, (10) 
where the operator D is defined by the relation (||) 
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Consider the T = case. One can use the low energy theorem for the derivative of the gluon 
condensate with respect to the quark mass 

(G^) = / d^x{G'iO)qq{x)) - ^^{qq) + 0{m,), (12) 



dniq 



where 0{mq) stands for the terms hnear in hght quark masses. Then one arrives at the foUowing 
relatiori0 

-{G^) + ]m = lm + ]m = m- as) 



Note that three fourths of the quark condensate stem from the giuon part of the nonperturbative 
vacuum energy density. Afong the same hues one arrives at the foUowing expression for the gluon 
condensate 

- De^a. = -^(4 - E + \ E ^M^)) = (G')- (14) 

g 9 q 

In order to get the dependence of the quark and gluon condensates upon T use is made of the 
Gell-Mann- Oakes-Renner (GMOR) relation (S \ {uu) \ = \ {dd)\) 

F^MI = + md){uu + dd) = [niu + md)T^ (15) 



Then we can find the following relations 

d T. d 



(16) 



The pressure of the massive relativistic pion gas has the form 

P^(T) = ~^ ^2^^1n(l_e--'(-)) = ^Vi / x^dxe- 
27r2 Jq 27r2 j-^nj^ 



Ax) 



n=l 

where uj-n{x) ~ -^/x^ + / and K2 is the Mackdonald function. Making use of (^,|l3|,|l6|) and 
( p^ ) one gets for the quark condensate 

p xHx 3M.r f.l M 

E 4^2j72y^ c^,(x)(e-'(-) - 1) 47r2F2^^n^^ T ^ 

In the chiral limit, M^r <^ T < Tc, the shift of the quark condensate is given by the standard 
temperature one-loop ChPT contribution to {qq)(T) 

^The equation (|l3[) contains corrections ~ mqd{qq)/dmq ~ 0{mq) which are negligible for u and d quarks. 
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Figure 2: The quark condensate S(T)/E as a function of the temperature. The sohd hue corre- 
sponds to Eq. (H^) . The dash-dotted hne corresponds to the three- loop ChPT formula (|l]) in chiral 
limit for Nf — 2. The dashed line is the numerical three- loop result of ChPT with non-zero quark 
mass ( M^=140 MeV ) from Ref.§. 



AS 



(21) 



E 47r2F2 ^ rfi 8F^ 

" n— 1 " 

In the opposite limit of low temperatures, T <^ Mt^, one obtains the following result 
AS 3AUT M , SM'J^T^^ 

E 4tt^F^ T ' 25/2^3/2^72^ ■ ^^^> 

The value of the quark condensate as a function of the temperature is depicted in Fig. 2. 
One can see that the gas formula (^0|) for the relativistic massive pious perfectly agrees with the 
numerical calculation of E(r) Q within the framework of ChPT with non-zero quark mass. At 
the temperature T = M-^ the deviation is only 2.5%. 

4. Now let us consider the temperature dependence of the gluon condensate within the frame- 
work of the present approach. Eqs. (M), (UJ) and (nsh yield 



(G^)(r) _ 24 M^T^ x^dx vIr- Ml 

(G2) " b (G^) Jo c^.(a;)(e-'(-) - 1) " M^^) n T 

For <t:T one gets 



(23) 



A(G2)_ 24M|T2^J_^_47r2M|T2 



(G2) b (G2) ^^n^ b (G2) 



In chiral limit, — 0, we have A(G2)(r) = 0, in agreement with the fact that a free gas of 
massless particles is conformally invariant. In the opposite nonrelativistic (low temperature) limit, 
T ^ Airr , one obtains 

A(G2) _ _24M|T _3^2^^M^^ 
(G2) " b (G2) T ^ b (G2) ■ 
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On the other hand, as stated earher, withm ChPT and for zero quark mass the dependence of the 
gluon condensate on the temperature arises only at the three-loop level due to the interaction of 
the massless pions @. In Refs. g the pressure of the Nf massless quarks has been calculated 
at the three-loop level with the result 



Nj=2 



ChPT 



30 ^ 42F2 



)2ln^} + Or°), 



(26) 



where the numerical value, relevant for the limit m„ — nid = at fixed is Ap ~ 275 MeV §. 
Then using ( [l0| ) , ( ITTl) and ([l6| ) one gets for the gluon condensate the result of Leutwyler 



{G'){T)chPT _^ 167r4 _ 1 

(G2) 1356F4(G2)^'' T 4^ 



(27) 



In Fig. 3 we show the temperature dependence of the gluon condensate at (G2) = 0.5 GeV*.^ 
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Figure 3: The gluon condensate {G'^){T) / {G"^) as a function of the temperature. The solid line 
corresponds to Eq.(|2|) at (G2) = 0.5 GeV*. The dash-dotted line corresponds to the three-loop 
result of ChPT in the chiral limit (Eq.(|3)). 

One can get an insight into the physical nature of small (G2) shift with the increase of the 
temperature. The smallness of this quantity is due to the large value of (G2) at T = as compared 
to the typical parameters of the thermal hadronic phase, A{G'^)/{G'^) oc M^T^ / {G'^) ~ 10^^ at 
A'U = 0.14 GeV, Tc = 0.17 GeV and (G2) = 0.5 GeV*. 

5. Within the above framework one can derive the thermodynamic relation for the quantum 
anomaly in the trace of the energy-momentum tensor of QCD. At low temperature the main con- 
tribution to the pressure comes from thermal excitations of massive pions. The general expression 
for the pressure reads 

Pn = TMMn/T), (28) 

3ln case = one has (Td/dT - 4)P = e - 3P = (e®^)T [i- 

^This value corresponds to the standard magnitude (^P^^) = 0.012 GeV* tedj. 
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where tp is a function of the ratio M^^/T and in case of the dilute gas it is given by equation ( p^ 
Then the following relation is valid 



d d 
(A - T Af^ 



2 " )P, = Ml^^^ 



Making use of (|g,|lO|), (|13|,22) and (gg) one gets 



dm„ 



where A{qq) = {qqjr — {qq) and A{G^) — (G'^)t — (G^)- In view of ( |i7| ) one can recast (| 
form 

Let us divide both sides of (|l|) by AT and take the limit AT 0. This yields 



d{G^] 
dT 



32^2 



djqq) 
dT ■ 



(29) 

(30) 
in the 
(31) 

(32) 



This can be rewritten as 



d{9s 



MA" 



(33) 



where (619^; 



= Y.'^qiqq) and 



9f , 
MM' 



dT dT 

a.^) /lQTia1){G'^) are correspondingly the quark and gluon 



contributions to the trace of the energy-momentum tensor. Note that when deriving this result use 
was made of the low energy GMOR relation, and therefore the thermodynamic relation ( ^2||3^ ) is 
valid in the light quark theory. Then at low temperature T when the excitations of massive hadrons 
and interactions of pions can be neglected, equation ( p3[) b ecomes a rigorous QCD theorem. This 
can be easily verified via direct calculation using (pO|), (p3|) and the GMOR relation. 

6. The present paper is devoted to the thermodynamic properties of QCD nonperturbative 
vacuum with two flavors at low temperature outside of the scope of perturbation theory. It was 
shown that at temperatures T < Af^ the relativistic massive pions can be treated within the dilute 
gas approximation. Analytic temperature dependence of the quark condensate perfectly agrees, in 
the low temperature region, T < Af^r, with the numerical calculations of {qq){T) obtained at the 
three-loop level of ChPT with non-zero quark mass 0. The gluon condensate slightly varies with 
the increase of the temperature, i.e. the situation is similar to the chiral perturbation theory. It 
was shown that the temperature derivatives of the anomalous and normal (quark massive term) 
contributions to the trace of the energy-momentum tensor in QCD with light quarks are equal to 
each other in the low temperature region. 

As it was mentioned above the pion plays an exceptional role in thermodynamics of QCD due to 
the fact that its mass is numerically close to the phase transition temperature while the masses of 
heavier hadrons are several times larger than Tc. This was the reason we did not consider the role 
of massive states in the low temperature phase. This question was discussed in detail in Ref. 0. 
It was shown there that at low temperatures, the contribution to (qq) generated by the massive 
states is very small, less than 5% if T below 100 MeV. At T = 150 MeV, this contribution is of 
the order of 10%. The influence of thermal excitations of massive hadrons on the properties of 
the gluon and quark condensates in the framework of the conformal-nonlinear a- model was also 
studied in detail in [ p3| . 
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